Fractal Substructure of a Nanopowder 
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The structural evolution of a nano-powder by repeated dispersion and settling can lead to char- 
acteristic fractal substructures. This is shown by numerical simulations of a two-dimensional model 
agglomerate of adhesive rigid particles. The agglomerate is cut into fragments of a characteristic 
size £, which then are settling under gravity. Repeating this procedure converges to a loosely packed 
structure, the properties of which are investigated: a) The final packing density is independent of 
the initialization, b) the short-range correlation function is independent of the fragment size, c) the 
structure is fractal up to the fragmentation scale £ with a fractal dimension close to 1.7, and d) the 
relaxation time increases linearly with t. 

PACS numbers: 45.70.-n,45.70.Qj,61.43.Gt,61.43.Hv 



The van-der-Waals attraction between nano-particles 
is much stronger than their weight. This is the reason 
why they agglomerate into ramified, often fractal struc- 
tures [lj. A well studied example are the agglomerates 
formed in a filter that collects nano-particles. Once the 
particle deposits have grown to micrometer size, they can 
be shaken off the filter fibers easily. Collecting these 
rather large aerosol flakes in a container leads to what 
is commonly called a nano-powder, a fragile assembly of 
partly sintered micrometer flakes made of nano-particles. 
Depending on the agglomeration process in the aerosol, 
as well as the influence of diffusion on the deposition pro- 
cess, the nano-powder will have fractal substructures Q. 
However it may be questioned, whether these are robust: 
Shaking, pouring, stirring, and all kinds of random treat- 
ments of the container will break the nano-powder up into 
fragments, presumably with a typical size determined by 
the prevailing shear forces and much larger than the pri- 
mary nano-particles, but not necessarily larger than the 
originally collected aerosol flakes. When allowed to set- 
tle, these fragments will reagglomerate, until the next 
perturbation breaks the nano-powder up again. 

Nano-powders can be enormously porous. Porosities of 
more than 90 % are common. Since many physical prop- 
erties, such as electrical conductivity, mechanical stabil- 
ity, or catalytic activity, are determined by the structure 
of the powder, it is important to know, whether there 
emerge robust generic structural features as a result of 
repeated fragmentation and reagglomeration processes. 
In this paper we present large scale simulation results 
for a simple two-dimensional model which shows such a 
development of a robust asymptotic structure. 

In our model the nano-particles are represented by up 
to 3 million discs with a narrow size distribution (10% 
variance). As initial state we take a densely packed ag- 
glomerate. Below we will show that the final structure 
is independent of the initial configuration. Then the fol- 
lowing procedure is repeated many times: First the ag- 



glomerate is cut with a square mesh into portions. The 
linear mesh size I can be viewed as the typical scale of the 
fragmentation process. A portion may consist of several 
disconnected fragments. These fragment flakes then set- 
tle as rigid bodies under gravity without taking adhesion 
forces with other particles into account. This is justified, 
if the flakes are sufficiently large, so that their weight ex- 
ceeds the van-der-Waals force between the nano-particles 
[Icj |. Brownian motion is neglected for the same reason. 
After this reassembly of the fragments the agglomerate 
is cut again with the square mesh, and so on, see Fig. [T] 

The only model parameter is the linear mesh size I. 
Note that the limit of small I {I k, particle diameter) cor- 
responds to non-cohesive primary particles. In any case, 
the cohesion forces are assumed to be weak compared 
to the fragment weight, but strong enough to assure the 
internal stability of the flakes. 

Algorithm: To obtain statistically significant results 
for the structure of nano-powders one has to consider 
systems with more than one million particles for many 
fragmentation-reagglomeration cycles. This is beyond 
the capability of Molecular Dynamics simulations. For 
our purpose we therefore generalized a model by Vis- 
scher and Bolsterli [3j originally intended for the sequen- 
tial deposition of macroscopic spherical particles (see 
0, IE S 0, Q for other applications): Each particle 
starts at a random position well above the already de- 
posited material (the configuration of which is regarded 
as frozen in) . Following gravity, it moves downwards until 
it touches the bottom of the container, where it sticks, 
or contacts another already deposited particle. In the 
latter case it moves, again following gravity, on the sur- 
face of the deposit until it either touches the bottom or 
finds a stable position in contact with the walls and/ or 
previously deposited particles (for more details see [9j). 

We generalized this algorithm in order to apply it to 
the fragments of a nano-powder. In each iteration step 
we inspect the portions which are cut out by the square 
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FIG. 1: Evolution of the packing of a nano-powder as de- 
scribed in the text, a) Initial packing generated by random 
sequential sedimentation Q|. The packing is cut by a square 
mesh into fragments (£ — 20); b) The fragments are consid- 
ered as rigid bodies and deposited (1 st generation). Again 
the packing is cut by a square mesh (here not shown); c) the 
fragments are deposited again (2 nd generation), and so on; d) 
3 rd generation; e) 4 th generation; f) 120*'' generation. 

lattice, Fig. [TJ with respect to their connectivity. A por- 
tion may decompose into several fragments (i.e. clusters 
of connected nano-particlcs). We deposit these fragments 
in a random sequence in the same way as described above. 
A fragment rolls down the surface of the deposit until the 
vertical projection of its center of mass falls in between 
two points of contact. As in the original algorithm, in- 
ertia is neglected, which in contrast to previous applica- 
tions is less of a problem here, because the dynamics of 
nano-particle flakes is usually strongly damped. In the 
following, lengths are given in units of the average parti- 
cle radius, masses in units of the particle mass, and time 
as number of fragmentation-reagglomeration cycles. 

Asymptotic filling height. The original Visscher- 
Bolsterly algorithm produces random dense packings of 
spheres without fractal substructures. Correspondingly, 
our generalization produces a packing of fragments that 
is homogeneous on scales larger than the fragmentation 
length t, as can be seen in Fig. [if. Surprisingly, however, 
the short range structure up to size i develops robust 
fractal properties. A first indication is given by the l- 
dependence of the filling height. 

Starting from a random dense packing of primary par- 
ticles the filling height increases towards a saturation 
value. Asymptotically, the powder adopts a very porous, 
statistically invariant structure, which is robust with re- 
spect to fragmentation at a fixed scale and subsequent 
gravitational settling of the fragments. Remarkably the 
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FIG. 2: Evolution of filling height starting from a random 
dense packing of height /io (indicated by the dashed line). The 
same asymptotic filling height is reached from above, if the 
particles initially form a single vertical needle (data marked 
by +) instead of a random dense packing. The full lines are 
fits according to Eq. (fl]). 
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FIG. 3: The asymptotic filling heigth h x grows as a power 
law hao(£) ~ £ a with mesh size I. The full line shows the 
best fit, a = 0.327. Inset: The relaxation time n c (l) increases 
linearly with mesh size. 

asymptotic filling height does not depend on the initial 
configuration: Starting with all particles arranged in a 
single vertical needle-like chain leads to the same value 
(see Fig. [2]). Except for the first point (the initial 
condition) the filling height h n at iteration step n can 
be fitted by an exponential approach of the asymptotic 
height, /too, with a relaxation time, n c , 

K = hoo(£) - (/looW - h )exp [-n/n c (t)\ . (1) 

The inset of Fig. [3] shows that 

n c (t)<xl z with z = l. (2) 

For the asymptotic filling height, a power law 

h 00 {t)<xl a with a = 0.327 (3) 

gives a very good fit(see Fig. [3|). 

This implies that the number of portions cut from the 
steady state configuration of a system of width L scales 
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FIG. 4: Pair correlation function (semi-log-plot) of the steady 
state structure. 
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FIG. 5: Average fragment mass as a function of the mesh size 
t in the asymptotic steady state. 



like N„ 



L/£ 2 oc £ a 2 . Consequently, the mass per 

d. 



portion is M/N p oc f' with 

dt = 2-a = 1.67 ±0.03. 



(4) 



One can interpret df as the fractal dimension of the 
structure on length scales smaller than the fragmentation 
length £, because the short range part of the pair corre- 
lation function g(r) does not depend on £ (see Fig. 
This means that the short range structure is indepen- 
dent of the overall density. For large distances, on the 
other hand, the pair correlation function approaches the 
overall density, which decreases with increasing I. 

On first glance, however, the fractal dimension Eq. (|3|) 
seems at odds with the fact, that the asymptotic aver- 
age fragment mass grows linearly with the mesh size £ 
(Fig. 03), suggesting, instead, that the fragments are effec- 
tively one-dimensional structures. This puzzle can only 
be resolved by assuming that the number of disconnected 
fragments Nf per mesh cell has itself a power law depen- 
dence on I; 
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As we are going to prove in a moment, the 
fragmentation-reagglomeration dynamics implies that 
the exponents (3 and df must be related by 

= df - 1. (6) 

Consequently the mass per fragment is linear in I: 

M £ d{ 

The proof of the scaling relation Eq. ([6|) is based on the 
steady state condition that the deposition of fragments 
reestablishes on average as many contacts as were cut in 
the preceding fragmentation step. An I x t- mesh cell 
contains oc £ di particles, of which oc £ dl ~ 1 are at the cell 
boundary. Hence, the number of particle contacts cut by 
the boundary of one cell scales as 
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FIG. 6: Normalized fragment mass distribution for different 
mesh sizes I. f(m) is the number of fragments of mass m 
divided by the total number of fragments for a given I. 



Each of the I® fragments forms two new contacts when 
deposited. Equating A cut = 2N{/N p gives Eq. ©. 

Fragment mass distribution. A detained understand- 
ing of the fragment properties is provided by the distri- 
bution of fragment masses, shown in Fig. [6] It reveals 
that that one must distinguish two types of fragments, 
large chunks at the upper end of the mass spectrum with 
a characteristic size m c , and scale invariant dust respon- 
sible for the power law part that is cut off by m c . Com- 
paring the mass distributions for different mesh sizes £ 
shows, that they can approximately be written in the 
form 



f(m,£)=m- T f 



m c {£) 



(9) 



(8) 



where the scaling function fix) is constant for x <C 1, 
goes through a maximum at x = 1, and has an approx- 
imately Gaussian tail for x 3> 1. The typical mass m c 
of the chunks has a power-law dependence on the mesh 
size, m c = 0.304 £ 1695 (Fig. [7|), the exponent being in 
good agreement with the value of df, Eq. ((4]). 

The evaluation of the dust exponent, r, is more diffi- 
cult, since the slope fitted to the power law part of the 
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FIG. 7: Chunk mass m c as a function of mesh size. Slope of 
straight line is 1.695. 




FIG. 8: Effective dust exponent r vs. l/£. 



mass distribution in the log- log-plot Fig. [6] decreases sig- 
nificantly with increasing I. An extrapolation of r for 
(r 1 / 2 — > 0, see Fig. El gives an estimate of r w 1.38. 
However, an independent method of determining t gives 
a larger value. It is based on the important observation 
that the width of the chunk-distribution is proportional 
to m c . Hence the fraction of chunks among the fragments 
vanishes like 

/chu, 1 ks«m^«^ (1 - T) , (10) 

because r is larger than 1. Hence, the normalization of 
/(m) for £ — ► oo implies that 

l=/(l)^m- T =/(l)C(r) (11) 

m 

with Ricmann's zeta-function at the argument r. Solving 
this equation numerically with /(l) ~ 0.36 (for the l- 
values we considered) gives an estimate r ~ 1.46. 

We have seen, that the overwhelming number of frag- 
ments are dust particles, apart from a vanishing fraction 
of chunks. However, these dust particles carry only a 
vanishing fraction of the total mass M. According to 



Eq. 0, 

N f 1 

m dust </(l)C(T-l)^oc- (12) 

vanishes for I — > oo. Essentially all the mass is in the few 
chunks. This explains, why the mass (essentially mass of 
chunks) per fragment (essentially per dust particle) has 
nothing to do with the fractal dimension. 

Now a consistent picture has formed: Each portion (or 
mesh cell) contains typically one chunk. The number of 
fragments per portion, which according to Eq. {5} scales 
like £@ , can thus be identified with 1 //chunks, which ac- 
cording to Eq. ([TO 7 ]) scales as £ < W r ~ 1 ). This shows that 
the fractal dimension of the chunks and the dust expo- 
nent are not independent of each other. Using Eq. ^ 
they obey the scaling relation 

d f {2 - t) = 1 . (13) 

For df = 1.695 this implies t — 1.41, in between the two 
r-values obtained above. 
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